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Status of weak matrix elements is reviewed. In particular, e'/e, B — > K tr y, Bk, Bb and Bb, are discussed 
and the overall situation with respect to the lattice effort and some of its phenomenological implications are 
summarized. For e'/e the need for the relevant matrix elements is stressed in view of the forthcoming improved 
experiments. For some of the operators, (e.g. 0$), even bounds on their matrix elements would be very helpful. 
On B —> K*j, a constant behavior of T2 appears disfavored although dependence of T2 could, of course, be milder 
than a simple pole. Improved data is badly needed to settle this important issue firmly, especially in view of 
its ramification for extractions of Vtd from B —> pj. On Bk, the preliminary result from JLQCD appears to 
contradict Sharpe et al. JLQCD data seems to fit very well to linear a dependence and leads to an appreciably 
lower value of Bk- Four studies of Bk in the "full" (71/ = 2) theory indicate very little quenching effects on Bk', 
the full theory value seems to be just a little less than the quenched result. Based on expectations from HQET, 
analysis of the B-parameter (Bhi) for the heavy-light mesons via Bhi = constant + constants'/m^ is suggested. 
A summary of an illustrative sample of hadron matrix elements is given and constraints on CKM parameters (e.g. 
Vtd/Vts), on the unitarity triangle and on x s /xd, emerging from the lattice calculations along with experimental 
results are briefly discussed. In quite a few cases, for the first time, some indication of quenching errors on weak 
matrix elements are now becoming available. 



Introduction 

The Lattice method is now in use for about a 
dozen years for calculating weak matrix elements 
[1,2]. Indeed there are by now about a dozen 
groups (APE, Columbia, FNAL, JLQCD, LANL, 
MILC, NRQCD, GKPS, UKQCD, Wuppertal. . .) 
involved in such efforts. In addition to these many 
groups there are also several individual attempts. 
It is abundantly clear that the hadronic matrix 
elements effort has evolved into a major activ- 
ity amongst the lattice community. This mas- 
sive amount of effort, though, is concentrated 
on a very few types of problems mainly on two 
and three-point functions. The problems that at 
first attracted many of us to the lattice, e.g. the 
AI =1/2 rule, requiring 4-point function calcu- 
lations, remain essentially unresolved to this day. 
Fortunately even for the simpler and limited class 
of problems entailing computations of 2- and 3- 
point functions we can make important impact 
to phenomenology. So the efforts are certainly 
worthwhile. 

One extremely attractive ("sexy") and rather 
unique feature of the weak matrix effort on the 



lattice, that has been recognized for a long time, 
is that it can have repercussions far beyond QCD. 
Since it is difficult to review all the work that is 
going on anyway and since the primary challenge 
facing the Particle Physics community at large, at 
present, is the quest for new physics, I will choose 
a few calculations below that tend, in my opin- 
ion, to enter into considerations that illustrate 
this feature. Here is the outline: 

1) e'/e; 2) B — > K*j and related matters; 3) 
B K ; 4) Bb', 5) Sample of Hadron Matrix El- 
ements Results; 6)Lattice + Experiment Con- 
straints on the SM; 7) Summary. 

For e'/e the theoretical underpinning as they 
pertain to the lattice have been greatly clarified 
in the last few years through the relentless efforts 
of Guido Martinelli and collaborators [3]. Unfor- 
tunately, there is acute lack of progress in the lat- 
tice calculation of the relevant matrix elements. 
One of the reasons for including this topic in my 
review is to stress this point and to emphasize 
that the experimental status implies an impor- 
tant window of opportunity for the lattice. 

Four groups [APE [4], UKQCD [5,6], LANL 
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[7], BHS [8]] have been working on the matrix 
elements for B — > K*j. In particular UKQCD 
has recently reported analysis of B — > K*j as 
well as of the decays B — > -kIvi [9] and B — > 
p£i<7 [10] that are all related through Heavy Quark 
Symmetry (HQS) [11,12]. 

On Bk new results are being reported using 
both staggered and Wilson fermions. Perhaps 
most interesting are those by the JLQCD group 
[13]. The Columbia [14] group is reporting a 
relatively high statistics study using dynamical 
fermions (i.e. nj = 2). Both of these works use 
staggered fermions. Preliminary results with Wil- 
son fermions from LANL [15], APE [16] and us 
[17] are also being reported. We have some pre- 
liminary results on Bk with nj = 2 "full QCD" 
configurations using Wilson fermions. 

Comparison of all the four unquenched results 
with their quenched counterparts shows that in 
each case B K (Q = Quenched) is just a bit bigger 
than B K } ~ 2 . 

We [17] are also presenting preliminary results 
for the B-parameters for B-mesons i.e. for bs and 
bd bound states including an exploratory study 
using n-f = 2, (3 = 5.7 unquenched configurations. 
Again we found that these parameters are a bit 
smaller (i.e. just a few per cent at most) in the full 
theory over their value in the quenched case. Our 
[17] data also shows clear evidence that Bi, s (= 
Bb s ) is a little (at most a few per cent) bigger 
than B bd (= B B )- 

Finally, I present a sample of the results for 
hadron matrix elements from the lattice that are 
of special importance to phenomenology. I then 
discuss the implications of these lattice results for 
weak matrix elements for constraining the param- 
eter space of the SM. For this purpose I briefly 
review the situation with regard to the Wolfen- 
stein parameters p, r\ as well as a, (3, i.e. two of 
the angles of the unitarity A. Numbers emerging 
for Vtd/Vt s as well as the implications for B s — B s 
oscillations are then briefly discussed. 

1. e'/e 

It is perhaps appropriate to begin the discus- 
sion with a brief reminder of the experimental 
status. After some quarter of a century of inten- 



sive efforts the experiments have produced two 
somewhat conflicting results: 

\e'/e\ = (23 ± 7) x 10" 4 NA31 (CERN) [18] 
\e'/e\ = (7 ±6) x 10" 7 E731 (FNAL) [19] (1) 

In addition to the disagreement, which is numer- 
ically quite significant, the interpretation of the 
results of the two experiments leads to strikingly 
different conclusions. E731 being consistent with 
zero means the Superweak Theory (SW) [20] can- 
not still be ruled out. On the other hand, if NA31 
is correct, not only the SW is clearly ruled out, 
it is also not inconceivable, as I will stress be- 
low, that the result is too high even for the SM 
and it presents an extremely important hint for 
physics beyond the SM! The choice between the 
two alternatives (SM or beyond) is the task of the 
theorist. It is in this context that lattice calcu- 
lations have an unusually significant role to play, 
as will be emphasized in the pages to follow. 

In passing it is also worth noting that improved 
experiments by the two groups are now well on 
their way on both sides of the Atlantic (NA48 at 
CERN and E832 (&TEV) at FNAL). These are 
promising to reduce the error significantly, i.e. 
from about .0007 at present down to .0001!, in 
the next few years. In particular, it is impor- 
tant to bear in mind that if NA31 is correct then 
both the experiments should be able to report a 
clear non- vanishing signal for e'/e even before full 
completion of the experiments i.e. even when the 
accuracy is only, say, .0002 or .0003 and not the 
aimed one of .0001. Thus the lattice weak matrix 
element effort may have a dramatic role to play 
on the fate of the SM in a few years time. 

Recall that e' is a measure of "direct" CP 
violation, i.e. occurring at the AS = 1 decay 
vertices as a difference between the amplitudes 
{■wk\H w \K°) and {itk\H w \K°). Thus 



(A •\H W \K L ) 
(Tr°Tr°\H w \Ks) 
(ir+ir-\H w \K L ) 
(ir+ir-\H w \K s ) 
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The calculation of the AS = 1 effective Hamilto- 
nian to the next-to-leading order (NLO) has been 
done by two groups [3,21]. The resulting expan- 
sion involves the usual ten four quark LL, and 
LR operators. 

In evaluation of the matrix elements of the left- 
right (LR) operators (i.e. 05_s) one needs the 
strange quark mass through, e.g. [3] 

(0 6 )vsa ~ {Ik - U)m A K /[m s {n) + m d (fi)] 2 (3) 

Consistent calculations of the matrix elements 
relevant to e'/e thus requires m s at a scale /x. It 
should be clear that the lattice approach has the 
significant advantage that such scale dependent 
calculations are feasible. 

Indeed, considerable progress has already been 
made in calculating m s (/x). In particular, Allton 
et al. [22] use input from many quenched simula- 
tions and to NLO determine: 



mf b (n = 2 GeV) = 128 ± 18MeV 



(4) 



In this calculation data from various quenched 
simulations with 6.0 < j3 < 6.4 with Wilson and 
Clover actions is used. From these data m s (a) as 
well as its tadpole improved counterpart rh s (a) 
are obtained through the relations: 



kc kr 



Ke 



These are then related to mf 5 (/i) to NLO accu- 
racy yielding (4). 

The error quoted in (4) is claimed to include 
systematics due to Clover versus Wilson, and 
due to boosted coupling [23] versus not boosted, 
m s (a) versus rh s a, and also scale uncertainties. 
While the error analysis in such a calculation is 
clearly rather difficult, the range of 128 ± 30 ap- 
pears quite safe. 

It is important to note that the scale, /x, 
most natural for typical lattice computations is 
~ 0(a _1 ) i.e. about 2-4 GeV. This is very for- 
tunate as the Wilson coefficients at such a rela- 
tively high scale can be systematically calculated 



by using improved perturbation theory. Indeed at 
/x<^ 1 GeV, which is more suitable for some other 
non-perturbative methods (e.g. Chiral perturba- 
tion theory (%PT)) Wilson coefficients of some of 
the important operators (e.g. O5 and 0e) are ex- 
tremely sensitive to the precise numerical value 
of (j, rendering such methods rather ill suited for 
reliable calculations [24]. 

Recall that, in the continuum, by using xPT, 
K — > 2tt matrix elements can be related to those 
for K — > 7r and K — > vacuum [25]. For lattice 
computations, this means that the much harder 
4-point functions can be obtained as linear com- 
binations of 3- and 2-point functions. Unfortu- 
nately, the method cannot be used with Wilson 
fermions as they do not respect chiral symmetry 
on the lattice [26]. It does work for the stag- 
gered case and in that approach it has become 
the method of choice [27]. 

They are several advantages of introducing "B- 
parameters" for these matrix elements, which 
measure the deviation from vacuum saturation. 
Thus: 

B i = (*\Q i \K)/(*\Q i \K) VS A (6) 

Several years ago, the staggered group made some 
progress in calculating B5 and Be [28]. They used 
P = 6.0 on 24 3 x 40 and 16 3 x 40 lattices. The 
results were: 



-Br , Br 



1 ± .1 



(7) 



These calculations, while pioneering, were rather 
crude. Comparison of the 16 3 and 24 3 lattices 
showed substantial finite size effects [27]; it is 
therefore not clear how large are volume errors on 
their bigger lattice although the general expecta- 
tion is that the 24 3 lattice should be sufficient. 
Perhaps a more serious problem is that the beta 
dependence of B5, Be as well as of just about 
every matrix element relevant to e'/e has never 
been studied. The only exception in this regard 
is Bk which enters e'/e indirectly (see below); for 
Bk the dependence on lattice spacing is quite sig- 
nificant especially for (3 < 6.0 [29]. So the num- 
bers above (7) must be regarded as rather ten- 
tative. Ironically, none of these matrix elements 
have received any attention in the last 5-6 years. 
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Figure 1. Theoretical expectations for e'/e [3,32], 
(solid 68% CL, dashed 95% CL). Experimental 
results [18,19] are also shown. 



With the improved computational resources that 
are now available many of these matrix elements 
should be calculable with appreciable precision at 
least in the staggered formulation. 

With Wilson fermions, in principle, the most 
attractive method is to attempt K — > 2tt "di- 
rectly" [30] (i.e. without reducing it to K — > 7r 
and K — > vacuum), wherein the use of the dis- 
crete symmetry CPS makes the theoretical un- 
derpinning very clean. I am quite hopeful that 
with today's computational resources this 4-point 
functions method may also work. Indeed, in this 
regard, the KEK method [31] for calculating mul- 
tipoint functions seems also to be a promising way 
to attempt. Finally, the use of the boosted cou- 
pling and tadpole improved perturbation theory 
[23] also tends to improve the chiral behavior of 
Wilson fermions. Thus the use of xPT to calcu- 
late the reduced 3- and 2-point functions also has 
a fair chance of working now. 

In their study Ciuchini et al. [3] use, based on 
existing, somewhat crude, lattice calculations and 
some "guesswork" (the ones with * on them): 



B K = .75 ± .15 B= 2 = 0-0.15 B* i = 1-6 

B B , 6 = 1.0±0.2 B\Ll* = 1.0 ± 0.2 B*£* = 1.0 ±0.2 

(note that hatted quantities are renormalization 
group invariant). Their key result is [32]: 

|e'/e| = (1.6 ± 1.9) x 10" 4 (8) 

How stable is this result to some of the inputs 
and B-parameters? To answer this the following 
variations were tried: 

a) Imposing an / B -cut of 200±40 MeV yields [32] 
\e'/e\ = (2.0 ± 2.3) x 10 -4 (9) 

b) Imposing a different /s-cut of 150 ± 25 MeV, 
motivated by the preliminary MILC results [33] 
you get [32]: 

|e'/e| = (1.9 ±2.1) x 10" 4 (10) 

c) Following the LANL group's preliminary re- 
sults [7,15], use ofi^ 2 = B 3 J 2 = 0.7 ±0. 15 yields 
[32]: 

|e'/e| = (1.7 ± 1.9) x 10" 4 (11) 

Thus one arrives at the conclusion that it is very 
difficult to get \e'/e\ > 10 x 10" 4 [3]. This, of 
course, leads to a contradiction (see Fig. 1) with 
NA31 which is the more "significant" of the two 
experiments, if you take the experimental errors 
at their face value. 

Here are some possible "loop-holes": 

1. 06 (or other matrix elements) suffer signif- 
icant final state interactions (FSI). In this 
regard, even bounds or clear demonstration 
of FSI would be helpful. 

2. m s (/x) is appreciably lower than the Allton 
et al. results [22]. This seems unlikely. 

3. NA31 could be "off". 

4. SM is "off". 

While the last conclusion is clearly the most dra- 
matic, it is entirely plausible. In this regard, it 
is important to note that this is SM with CP vi- 
olation. As is well known the likelihood of the 
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failure of the SM to account for CP violation is 
not negligible. After all theoretical investigations 
tend to strongly suggest that SM cannot account 
for baryogenesis [34]. 

It is clear then that there is a significant win- 
dow of opportunity for the lattice in conjunction 
with new experiments now in progress that will 
start to produce improved results in the next 2-3 
years, to decide the fate of the SM. 

2. B -> K*j and Related Matters 

Recall that radiative B decays, exclusive (e.g. 
B — > K*y, fry. . . ) or inclusive (i.e. b — > sj), are 
important tests of the SM and they are very sen- 
sitive to new physics [35]. Also the dynamics of 
the exclusive modes (i.e. B — > p(K*) + 7) is im- 
portant to understand as it controls extraction of 
Vtd from these decays. For this purpose a precise 
understanding of SU(3) breaking in these decays 
is also very important. 

We recall the CLEO measurements [36]: 

BR(B K*j) = (4.0 ± 1.7 ± 0.8) x 10" 5 

BR(b -> sj) = (2.32 ± .57 ± .35) x 10 -4 (12) 

From the ratio of these two BR we get [Rk* = 
BR{B K*j)/BR(b sj)} 

Rk? ~ 20 ± 10% (13) 

While these experimental observations are ex- 
tremely important, their quantitative improve- 
ments are needed for a meaningful confrontation 
with theory. 

Recall that the hadronization ratio, Rk* , is the 
quantity of direct importance to the lattice. In 
this regard, the experimental signal for B — > jK* 
need not all have a short distance origin [37,38]. 
The first point is that the lattice calculations of 
B — > jK* is, by construction, the matrix element 
of a very well defined (short-distance), two-quark 
operator, O7 [8,39]: 

O r = -em i s a a!? +A b ol F> lv /AS* 2 (14) 

Furthermore, in identifying the matrix element 
of O7 : (B\Or\jK*) with the experimental obser- 
vation of B — > jK* , one is also making a tacit 



assumption that the spectator model works. In 
other words one is assuming that the light (spec- 
tator) quark in the B-meson is not playing any 
important role in B — > jK* decay. 

There are compelling phenomenological rea- 
sons due to V u i,V us << V-t s for believing that the 
second assumption above is very safe. (Indeed 
this is not the case for the related decay B — > fry, 
especially for B~ — > p~j [37,38]). 

However, the first assumption should not be 
taken for granted. The point is that experimen- 
tally it is known that [40] 

BR(B^iPK*) = (1.58 ± .28) x 10" 3 
BR(B — > ip' K*) = (1.4 ± .9) x 10" 3 (15) 

which are both substantially larger than B — > 
jK* . Quantum mechanics tells us that some 
fraction of the off-shell will contribute to 

B — > jK* . Such a contamination is an example 
of "long-distance" contribution. It is extremely 
difficult to make reliable estimates of such LD 
contributions. Since the experimental result for 
B — > jK* contains all contributions — short and 
long-distance — whereas the lattice calculation, by 
construction, is SD — the difference between the 
two is a measure of the LD contributions. 

There are by now 4 lattice groups [4-8] 
(UKQCD, APE, LANL, BHS) who have been 
studying B -> jK* on the lattice. At LAT'94, in 
his review, Guido Martinelli [1] emphasized that 
the crucial issue was q 2 dependence of the relevant 
form factors (T2 or Ti). The point is that if T2 
tends to be constant with q 2 then Rk* tends to be 
about 25-35% whereas if T2 is pole-like then Rk* 
is appreciably smaller, i.e. about 5-10%. There- 
fore, it is clearly very important to resolve this 
issue. 

On the lattice separate calculation of the ma- 
trix element of the vector and the axial piece of 
the operator is a better strategy as then one has 
at hand one very good check, namely: 

T 2 (0) = Ti(0) (16) 

Thus 

<#*(!,, AOIigflfr)) = 2e llvXa r 1 v (k) P x k°T 1 {q 2 ) 
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(K*( V ,k)\A,\B(p)) 



v ■ q 



V 

(17) 



We recall that at the end point (i.e. q 2 = g^ ax = 
(rrihi — my) 2 ) T± and T3 do not contribute. Be- 
sides at q 2 = g^ ax both the initial and final 
mesons are at rest requiring no momentum injec- 
tion rendering the lattice calculation of T2(g£ lax ) 
very clean. Of course, for q 2 ^ g^ ax momentum 
injection is required. Also although q 2 = 0, the 
point of direct physical interest, is not accessible, 
current simulation parameters do allow q 2 /m 2 ll to 
be quite small i.e. g 2 /m^<^0.1. 

2.1. Implications of HQS[4,5,8,9,11,12] 

HQS provides very useful scaling laws. Thus: 



T 2(g max ) ~ 
In addition 

„2 



1/2 



m 



hi 



l hl 



1 



fl8) 



fl9) 



for mv/mht — > 0. Thus if we parameterize the 
g 2 -dependence of form factors as 



T lt2 (q 2 ) = T lt2 (0)/(l-q 2 /m 2 u r 



(20) 



We see that HQS scaling laws suggest that if T 2 
tends to be constant (i.e. n = 0) then Ti must 
behave as a simple pole (i.e. n = 1). On the 
other hand, if T 2 (q 2 ) falls as a simple pole (i.e. 
n = 1) then Ti should fall as dipole i.e. n = 2. 

It must be stressed though that these expecta- 
tions based on HQS are valid only for my /mht — > 
0. In particular, if the lattice heavy-light mass ~ 
2 GeV (especially when the simulations are at 
(3 = 6.0) then mx* /mht ~ -5 which may not be 
small enough for HQS to have set in. 

With regard to constant versus a pole behavior 
of T 2 (q 2 ) two comments are in order. First it 
is difficult to justify a constant behavior of the 
form factor as one is dealing with a bound state 
to bound state transition. Of course, in principle, 



an approximately constant behavior could result 
from the cancellation of the amplitudes over two 
or more resonances. The second point to bear in 
mind is that it is easier to distinguish between 
a pole versus a constant as the heavy-light mass 
gets heavier. If you consider the ratio: 



nqL^)/T(Q) = [l-q 2 ma Jm 2 hl ]- 1 



(21) 



then for m^i = 1.5 GeV, ?max/ m w ~ -1 an d 
T (9max)/ T (°) = i- 11 - Tnis means at m hl ~ 1.5 
GeV to distinguish a constant from a pole re- 
quires the data to have an accuracy a lot bet- 
ter that 11%. On the other hand if mht is 
heavier, say 4.5 GeV then ?max/ m w ~ -6 an d 
T (9max)/ T (°) - 2 - 7 wnicn is appreciably differ- 
ent from 1 so that even with a 30% accuracy one 
may be able to decide between a constant versus 
a pole behavior. 

Table 1 lists some of the characteristics of the 
four groups and their findings with regard to q 2 
dependence of form factors. Curiously LANL and 
APE using (3 = 6.0 tend to favor a constant be- 
havior and UKQCD (/3 = 6.2) and BHS (/3 = 6.3 
and 6.0) favor a pole-like behavior for T 2 . Perhaps 
this is related to the comments in the preceding 
paragraph as clearly the heavy-light masses ac- 
cessible at (3 = 6.0 are less than at (3 = 6.2 and 
6.3. 

Perhaps the best strategy for deciding between 
the two extreme options, i.e. T 2 constant versus 
T 2 pole-like, is to examine the ratio 



^(sLJ/T^O) 



(22) 



(Recall that T 2 (0) = Ti(0)). Using the j3 = 6.0 
and 6.3 data of Ref. 8 [extracted from Table III 
of that ref.] one finds that iZ C onst tends to be 
larger than 1 especially for the (3 = 6.3 data point. 
Indeed fitting the data to a constant one finds 
-Rconst = 1.39±.14 i.e. over 2cr away from 1. Next 
we examine the same data slightly differently, as: 



R 



pole 



"hi 



r 2 (9Lx)m(o) 



(23) 



Now numerically one finds, iZ po i e = -98 ± .10. 
Thus at least this data set tends to favor the sim- 
ple pole behavior for T 2 and tends to disfavor a 
constant behavior fairly strongly. 
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Table 1 
For B - 



K*y, the groups, their lattice parameters and results. 









#of 




T 2 




Rk-% 


Group 




Size 


Configs 


Action 


tends to favor 


T2-const 


T 2 -pole 


LANL [7,15] 


6.0 


32 3 x 64 


100 


Wilson 


const 


27 ±3 


4-5 


APE [4] 


6.0 


18 3 x 64 


170 


Clover 


const 


31 ± 12 


5 ± 2 


UKQCD [5] 


6.2 


24 3 x 48 


60 


Clover 


pole 


35^ 


13+ 14 


BHS [8] 


6.3 


24 3 x 61 


20 \ 












6.0 


24 3 x 39 


16 


Wilson 


pole 


N/A 


6.0 ± 1.2 ± 3.4 




6.0 


16 3 x 39 


19 J 











Recall that HQS tends to relate B 
factor to B — > tt [5,9,11,12]. 



nB — K* 



B->K* 



ft 



K*j form 



(24) 



Similarly form factors for B — > K*j are related 
to the axial current form factors for B — > plv 
[5,9,11,12]: 



2T; 



B->K* 



(25) 



where SU(3) is assumed. 

Recently UKQCD has studied form factors for 
B -> irlv [5] and for B -> plv [9]. Their data 
for B — > ttIv quite strongly favors a pole (dipole) 
dependence for /o(/i) over a constant (pole) one. 
Similarly their data for B — > plv favors a pole 
behavior over a constant one for A\. Finally 
FNAL group [41] also has a preliminary result 
for the /+, f° form factors for B -> irlv (f3 = 5.9, 
16 3 x 32, 100 configs.). Their data has very good 
statistics and in particular a constant behavior 
for fo is strongly disfavored. 

2.2. More HQS tests on form factors 

UKQCD [9] has studied the HQS expectations 
on ratios of form factors. They find that HQS re- 
lation [^4i/2T2 = 1] holds to about 20% accuracy 
even at mj. However, V/Ai = 1 fails at rap and 
works quite poorly even at mj (at 30% level). 

BHS [8] studied the HQS relation 



Ci + C 2 /m u 



(26) 



Their data in the range of 2<mu <^4 GeV gave 
a reasonably good fit to this relation. 



2.3. Summary on B — > K*j and related 
matters 

1. A constant behavior for T 2 (or for Ai) ap- 
pears to be disfavored. 

2. Of course a dependence milder than implied 
by the simple pole ansatz for T 2 cannot at 
this point be ruled out. More data is needed 
to examine this possibility critically. 

3. Various tests of HQS seem to work quite 
well for form factors for mht ^ tid- A no- 
table exception is V/Ai. 

4. Finally a suggestion for analysis. Since 
HQS makes a very useful prediction for 3- 
point function e.g. for T 2 : 

y 'mhiT 2 (q' 'max) = Consti+Const2/?n/a+- • 

as well as a similar relation for the two point 
function i.e. the decay constant fht- This 
suggests that it may be better to analyze 
the data as ratios (3-point functions over 2- 
point functions) so that 72(g^j ax ) / fhi could 
have improved convergence with respect to 
rriht- The additional advantage would be 
that the lattice should be able to provide 
more accurate predictions for such ratios. 

3. B K 

We recall that Bk defined as: 



B K = (K\[s 7li (l - 75 )d] 2 \K) 

/(K\[s lf ,(l- l5 )d] 2 \K)\ VSA 



(28) 
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controls the CP violation in the neutral kaon com- 
plex through the parameter e. As is well known 
this parameter (e) has the unique feature that it 
is the one and only manifestation of CP violation 
seen so far. What we can learn from the mea- 
sured value of e about the origin of CP violation 
is limited by how well we know Bk ■ 

3.1. Bk with Staggered Fermions 

At LAT'93 Sharpe [42] reported significant 
progress in pinning down Bk with staggered 
fermions: 



B K DR (2 GeV) = .616 ± .020 ± .017 



(29) 



This extraction of Bk was based on Sharpe's 
proof that the lattice spacing errors occur only 
at 0(a 2 ) and not at 0(a). In view of the signifi- 
cance of Bk it is clearly important to continue to 
scrutinize it with both staggered and with Wilson 
fermions. 

I will report briefly on the following results: 

1. JLQCD [13] (staggered, quenched). 

2. Columbia group [14] (staggered, "full" 
QCD i.e. n s = 2). 

3. APE [16] (Clover action). 

4. LANL [15,7] (Wilson). 

5. Bernard and Soni [17] (Wilson, see also 
Lat'94). 

JLQCD group is reporting perhaps the most in- 
teresting result on weak matrix elements at this 
year's symposium. They have data with reason- 
ably high statistics at four values of (3. Their 
simulation parameters are listed below: 



5.85 
5.93 
6.0 
6.2 



Size 
16 3 x 32 
20 3 x 40 
24 3 x 40 
32 3 x 64 



# of Configs. 
60 
50 
50 
40 



a 



-l 



GeV 
1.34(3) 
1.58(4) 
1.88(5) 
2.62(9) 



0.80 
0.75 
0.70 
0.65 
0.60 
0.55 
0.50 
0.45 




O non-invariant 
• invariant 



0.3 0.4 
ma 



Figure 2. JLQCD fit of B K [NDR, 2GeV] to a 
linear dependence on the lattice spacing, i.e. m p a 
for both gauge invariant and non-invariant oper- 
ators; see also [13]. 



data excludes the form A+Ba 2 \ If true this result 
would contradict Sharpe's proof [42]. 

Fig. (3) shows the JLQCD data along with 
Sharpe et al. (marked "GKPS") at f3 = 6.0, 6.2 
and 6.4. The GKPS data tends to lie systemati- 
cally a little bit above the JLQCD one, although, 
a direct comparison of the two data sets (i.e. with- 
out renormalization factor etc.) shows that the 
discrepancy is quite reduced [13]. 



0.45 















• 


















O JLQCD 








• GKPS 



0.0 



0.1 



0.2 



0.3 0.4 
ma 



0.5 



0.6 



Fig. (2) shows their key result. Their data for Bk 
(2 GeV) fits very nicely linearly with a (the lattice 
spacing) i.e. to the form A + Ba. In contrast their 



Figure 3. Comparison of JLQCD data for 
gauge non-invariant operators with Gupta et al. 
(GKPS) [42]; see also [13] 
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The most striking point of the JLQCD result 
[13] is that their data yields such a very high % 2 
(~ 25) to the fit of the form A + Ba 2 . If their 
present result withstands further scrutiny it will 
be a clear contradiction to Sharpe's proof. How- 
ever, two obvious items that need to be checked 
before drawing strong conclusions are: 

1. Fits of the form A + Ba 2 + Ca A . 

2. Perhaps one should discount data at 5.7, 
5.9. Then JLQCD would need to run at 
additional /3's which they are planning to 
do anyway. So these questions will be ad- 
dressed in the next few months. Meantime 
their preliminary result is [13]: 

B K (2 GeV) = .497 ± .008(stat) 

±.014(0( ff 4 )) (30) 

This central value of Bk is significantly lower 
than that of Sharpe et al. If it is proven correct 
then it will clearly have very important implica- 
tions for phenomenology. Also the JLQCD value 
of Bk tends to be somewhat lower than implied 
by large N [43] — though not in disagreement (due 
to the size of the error (~ .05) that is usually 
quoted in the large N result) and would, on the 
other hand, be a little closer to the lowest order 
chiral perturbation theory value [44]. 

The Columbia group [14] is reporting a very 
interesting measurement of Bk in "full" QCD i.e. 
with r} } = 2, am = 0.01, (3 = 5.7 on a 16 3 x 32 
lattice with 250 configurations. Their key result 
is: 

B£ DR (/x = ir/a) = .659 ± .063 (31) 

This result can be compared with the quenched 
(staggered) one at (3 = 6.0 

1. Sharpe et al. [42] give B K DR (2 GeV) = 
.707 ± .008. 

2. JLQCD gives B K DR (2 GeV) ~ .69 ± .01 
(this is read off of their figure). 



The Columbia study is the third attempt to 
examine Bk with dynamical quarks using stag- 
gered fermions. The previous two were the work 
of Kilcup [45] and that of Ishizuka et al. [46]. 
In all three (3 = 5.7, nj = 2, with the lightest 
dynamical quark mass am = 0.01, is used. All 
three studies indicate very little difference with 
the quenched result at (3 = 6.0. So the general 
expectation that most of the difference between 
full QCD and QQCD can be accounted for by a 
shift in the lattice spacing seems to hold to a very 
good approximation. We will come back to this 
point a littler later on. 

3.2. Bk with Wilson Fermions 

As is well known, with Wilson fermions, cal- 
culation of Bk becomes rather problematic due 
to the fact that Wilson fermions do not respect 
chiral symmetry [26, 47]. Three approaches are 
currently being used to address to this problem: 
1) The LANL group is calculating the 4-quark 
matrix elements for several values of momentum 
transfer [7,15]. The non-chiral part can then 
be subtracted non-perturbatively. 2) The APE 
group [16] is using the Clover action and also 
evaluating the weak coupling corrections using 
a very interesting non-perturbative method. As 
Fig. (4) shows this seems to alleviate the prob- 
lem at least to some extent. 3) We [17,48] use a 
method wherein the mixing coefficients amongst 
the 4-quark operators are allowed to vary to sys- 
tematically restore chiral symmetry. 

The LANL [15] group is calculating the ma- 
trix elements on their 32 3 x64 f3 = 6.0 lattice 
(150 configurations) at 5 values of momentum- 
transfer: p = (0, 0, 0); (1,0,0); (1,1,0); (1,1,1) and 
(2,0,0). To 0(p 4 ) xPT predicts [7,15]: 

(K o {p 2 \0 As=2 \K o { Pl )) = a + (3m 2 K + 7Pi ■ p 2 

+<5im 4 f + 6 2 m K Pi ■ P2 
+<5 3 (pi-p 2 ) 2 (32) 

with pi = [rriK , 0, 0, 0). Here a, (3, 8i are pure 
lattice artifacts due to mixing with the wrong 
chirality operators. The coefficients 7, 82 and 63 
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Figure 4. Chiral behavior of the matrix ele- 
ments of the AS = 2 operator + as a func- 
tion of X = 8/3f K M K /\(0\sj 5 d\K)\ 2 . Results 
obtained by renormalisation in "standard" pertu- 
bation theory (SPT), in "boosted" perturbation 
theory (BPT) and non-pertubatively (NPT) are 
compared; see also [16] 



contain lattice artifacts [7,15]: 



B 



K 



8/i 



[7 + fcmjr + hm K (Ei + Ej) 



(33) 



where E' s are the energies corresponding to dif- 
ferent momenta. After subtractions they find at 
f3 = 6.0 [15]: 



B£ DR (2.3 GeV) = .74 ± .10 



(34) 



which is in very good agreement with the stag- 
gered result [42]: .707 ± .008. 

In our analysis the data set in Table 2 is being 
used, where 5.7jr stands for rjp = 2, am = 0.01 
lattices borrowed from the Columbia group [49]. 
Comparison of (3 = 5.7 p with (3 = 6.0 (quenched) 
16 3 x 39 lattices shows that the errors on Bk 
due to quenching are small and are within one 
sigma of the statistical and systematic errors due 
to other sources. We will take the difference in 
the two central values (i.e. 0.05) as estimate of 
quenching errors. 

A linear fit (i.e. of the form A + Ba) to our data 
yields: 

B K =0 {2 GeV) = .58 ± .06 ± .01 ± .05 (35) 



where the errors are due to 1) statistics and sys- 
tematics of chiral symmetry restoration, 2) finite 
size effects and 3) quenching. 

3.3. Summary of Quenching errors on Bk 

It is very interesting to note that all four [see 
Table 3] studies [14,17,45,46] of the effects of 
quenching (three with staggered and one with 
Wilson fermions) indicate that effects of fermion 
loops are small (i.e. 10%) and within sta- 
tistical and systematic errors. However, curi- 
ously enough, in all 4 cases the "full" theory re- 
sult seems to be systematically a little below the 
quenched result. Thus, in each case the ratio 



Rfl > 1, where Rji 



(2 GeV) 



This may 



(2 GeV) 

mean that fermion loop reduce Bk by a few %. 

4. Bb- B-parameter for Heavy-Light 
Mesons 

On this topic I will mainly be reporting our own 
results [17]. I want to take the opportunity to re- 
iterate that our first study of the B-parameter for 
heavy lights was indeed quite correct [50]. Work- 
ing at (3 = 5.7-6.1 although we had used rather 
heavy quark masses (am > 1) without using the 
Kro-Mac norm [52], this does not affect the B- 
parameters as they are ratios in which that norm 
just cancels. We recall the result of that study: 



B bd (2 GeV) = 1.01 ± .15 



(36) 



i.e. that for the heavy-light mesons vacuum sat- 
uration works to a very good approximation 
[50,51,53,54]. Indeed since the reduced mass of 
the heavy-light meson is of the same order as the 
light meson (e.g. kaon) the agreement with vac- 
uum saturation is non-trivial and was not antic- 
ipated. In any case the precise value of the B- 
parameters for the two B-mesons (Bd and B s )aie 
extremely important as they enter in constraining 
the SM CKM parameters through the unitarity A 
(see the following section). 

Analyzing our data we find [17] that Bh s (light 
quark held fixed near the strange quark mass) is 
increasing with rrihs and for rrihs > 3 GeV, Bh s 
is within about 10% of unity i.e. the vacuum sat- 
urated value. Note in particular that the data 
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Table 2 

Bk with Wilson fermions using the method Ref. 47 [see also Ref. 17] 



P Size Configs a^/GeY B K (fi = 2 GeV) 

5.7 16 3 x 33 60 1.17 .81 ± .03 

6.0 16 3 x 39 60 2.2 .66 ± .08 

6.0 24 3 x 39 40 2.2 .67 ±.07 

6.3 24 3 x 61 30 3.21 .59 ± .07 

6.5 32 3 x 75 40 3.9 .70 ± .05 

5.7 F 16 3 x 39 49 2.28 .61 ± .07 



Table 3 



Effects of Quenching on Bk 



Study 


_iC' =u (2GeV) 
R ' 1 -bV~ { 2 GeV) 


Remarks 


Kilcup [45] 


~ 1.09 ± .07 


Staggered 


Ishizuka et al. [46] 


~ 1.03 ± .05 


Staggered 


Columbia Group [14] 


~ 1.06 ± .10 


Staggered 


Bernard & Soni [17] 


~ 1.08 ± .15 


Wilson 



set [see Table 1] contains also the unquenched 
("full" QCD) configurations borrowed from the 
Columbia group [49]. This is the first attempt 
at studying the difference between the quenched 
and the unquenched lattices for By. 

Since gluon interactions with heavy quarks can 
be expanded in powers of 1/ my heavy quark ef- 
fective field theory suggests that the Bht may be 
expanded in inverse powers of my : 

B P U = Co + Cf /my + • • • (37) 

where the superscript P is included anticipating 
the /3-dependence on a lattice calculation. From 
the lattice perspective, analysis of the data, along 
such an expansion, should be very useful in at- 
taining precision. 

Fig. (5) shows Bhs [2 GeV] versus 1/my from 
our data [17]. We fit for Cq and cf for each 
lattice and typically we find 

~ 1.04 ±.05 (38) 
Cf ~ -(.4±.l)GeV (39) 

In particular, the negative sign of cf is unam- 
biguous and furthermore its numerical value of 
about 400 MeV is characteristic of the mass scale 
of the "brown muck" [11]. Thus you see that 



at the B-meson mass the C\ term would hardly 
make a correction of ~ 10% on the asymptotic 
value of the B u [2 GeV]. 

Using such a parameterization we have calcu- 
lated B for the B s meson (Mj s = 5.3 GeV) and 
we find (see Fig. 3) [17] 

B bs (2 GeV) = .97 ± .05(statistical) ± .01 (finite size) 
±.02(scale - breaking) ± .02(fitting) 
±.04(quenching) (40) 

Note again that B hs in "full" QCD (i.e. n } = 2, 
P = 5.7, am = 0.01) is a bit lower (< .05) then 
its quenched counterpart at P = 6.0. 

Finally the SU(3) breaking in the B-parameters 
is very important (see below) for relating B s -B s 
and Bd-Bd oscillations as well as for extracting 
V-td from the ratio of the oscillation parameters for 
the two mesons [55]. For such phenomenological 
applications we have studied the SU(3) breaking 
ratio: 

Rsd = Bts/Btd (41) 

shown in Fig. (6). The following qualitative fea- 
tures are seen: 

1. Bi s > Bid 
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1.0 



0.8 — 



0.6 1 — 
0.0 



1 hs 



o -p 



#3 



o ,3=6.0 16 J x39 
+ ,3 = 5. 7F 16 
□ ,3 = 6.0 24 3 x39 
o ,3 = 6.3 24 3 x61 



39" 



0.2 



0.4 



0.6 



Figure 5. B-parameter for heavy-light (light be- 
ing roughly at the strange-quark mass) i.e.Bhs 
vs. 1/mhs', see also [17]. Vertical line indicates 
the location of the B s meson. 



2. R„ 



1 as mhs — > co 



3. > R? d 



Numerically we find 



R n s f 2 

sa 



R sd = 1.01 ± .02(st) ± .02(sy) ± .03(qu) 

1.02 ±.04 (42) 



sa 



where errors due to statistics, systematics and 
quenching are as indicated. 

5. Sample of Hadron Matrix Elements Re- 
sults from the Lattice. 

Table 4 shows an illustrative sample of results 
for hadron matrix elements obtained from the lat- 
tice (recall that Ba = Bb, -B&s = Bb,', these 
notations are being used interchangeably). Few 
remarks are in order: 

1. On the B paramters since the lattice 
method is evolving into a precise method for 
calculating these quantities at a scale > 2 
GeV we are opting to quote their values at 




0.9 



o ,3 = 6.0 16 J x39 - 
+ ,3 = 5. 7F 16 3 x39 
□ ,3 = 6.0 24 3 x39 - 
Q | ^ = 6.3 24 3 x61 

4 5 



m hs (GeV) 

Figure 6. The SU(3) breaking ratio, R s d (i.e. 
Bbs/Bbd) vs. m hs ; see also [17] 



2 GeV. The renormalization group invari- 
ant counterpart (often denoted by Bk) can 
be obtained from Bk (2 GeV) by multiply- 
ing with the Wilson coefficient to a specified 
order [56]. 

2. The 90% CL summary that is given for Bk 
still uses the central value of Sharpe et al. 
[42] due to the preliminary nature of the 
JLQCD [13] result; the errors on Bk are 
increased over those given in previous such 
summaries [55] to reflect the discrepancy 
between the two. 

3. Errors due to quenching are indicated, in 
the Table, wherever applicable, by the letter 
(Q) after the error. 

6. The Noose: Lattice + Experiment Con- 
straints on the SM [55,65,66]. 

An important mission of the weak matrix el- 
ement effort on the lattice, in conjunction with 
experimental results, is to deduce reliable con- 
straints on the SM. In particular, much attention 
is focused on the Wolfenstein parameters p, r\ or 
equivalently the angles a and (3 of the unitarity 
triangle [67,40]. From Table 4 using the lattice 
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Table 4 



Illustrative Sample of Hadron Matrix Elements from Lattice QCD 



QUANTITY 


VALUE 


AUTHORS (REMARKS) 


B K (2 Gev) 


.616 ± .020 ± .017 


Gupta, Kilcup, Sharpe (Staggered) [42] 




.65 ± .15 


ELC (Wilson) [58] 




.58 ± .06 ± .05(Q) 


Bernard, Soni (Wilson) [48,17]. 




.50 ± .01 ± .05 


JLQCD [13,57] 


B K (2 GeV) 


.62 ± .15 


Most likely = 90% CL [59] 


Btd (2 GeV) 


1.0 ± .15 


Bernard et al. [50] 




1.16 ± .07 


ELC [58] 




.96 ± .06 ± .04(Q) 


B + S [17] 


Bbd 


1.0 ± .15 


Most Likely (90% CL) [59] 


Bbs/ Btd 


1.01 ± .04 


B + S [17] 




1.08 ± .03 ± .08 


Bernard, Labrenz, Soni [60] 


fa/ (MeV) 


174 ±26 ±46 


Bernard, et al. [50] 




190 ± 33 


Degrand, Loft [61] 




210 ±40 


ELC [58] 




185+4+_ 4 7 2 


UKQCD [62] 




208 ± 9 ± 32 


BLS [60] 




182 ± 3 ± 9 ± 22(Q) 


MILC [33] 


fa/ (MeV) 


197 ± 25 


Most Likely (90% CL) [59] 


faj (MeV) 


222 ± 16 


Degrand, Loft [61] 




234 ± 46 ± 55 


Bernard, et al. [50] 




230 ± 50 


TPT C Tc;fil 




212 ±4+_ 4 . 6 


TTTt*"or , 'n Tfioi 
ujw^^jj L J 




230 ± 7 ± 35 


riLb [bOJ 




198 ± 5 ± 10 ± 19(Q) 


MTT c r^^i 


/ Ds /MeV 


221 ± 30 


Most Likely (90% CL) [59] 


fa/ (MeV) 


205 ±40 


ELC [58] 




160 ± 6±H 


UKQCD [62] 




187 ± 10 ± 37 


BLS [60] 




151 ± 5 ± 16 ± 26(Q) 


MILC [33] 


fa/ (MeV) 


173 ±40 


Most Likely (90% CL) [59] 


fa s l (MeV) 


194+jtf 


UKQCD [62] 




207 ± 9 ± 40 


BLS [60] 




169±7±14±29(Q) 


MILC [33] 


f B J (MeV) 


201 ±40 


Most Likely (90% CL) [59] 


/b s / /b 


1 22+ 04 


UKQCD [62] 




1.11 ± .02 ± .05 


BLS [60] 




1.22 ± .04 ± .02 


FSG [63] 




1.10 ± .02 ± .04 ± .08(Q) 


MILC [33] 


/b s / /b 


1.16 ± .10 


Most Likely (90% CL) [59] 


fix. = r( r B ^;' = 6.0 ± 1.2 ± 3.4% 


BHS [8] 




13li 4 % 


UKQCD [5,9,64] 




5 ±2% 


APE [4,64] 




4-5% 


LANL [7,15,64] 
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value of Bk with the kaon CP violation parame- 
ter e leads to one set of curves (see Fig. (7)) which 
enclose the allowed domain of the SM. Similarly, 
the B parameter for B-mesons and fs with the 
experimentally measured B-B mixing parameter 
Xd, yields two more curves enclosing the allowed 
area for the SM by this set of considerations. Fi- 
nally, the experimental observations [40] of semi- 
leptonic B decays and the B-lifetime yield y*^ 
(we use V ub /V cb = .08 ± .02 and V cb = .04 ± .005) 
leading to the two concentric circles on Fig. (7). 




Figure 7. Constraints on the p— r\ plane. Concen- 
tric circles (solid) are due to V ub /V cb ; e and Bk 
yield dashed curves and Xd with /b-Bb determine 
the dotted domains. The remaining allowed re- 
gion (90%CL) is indicated (solid). 




ALPHA 

Figure 8. The remaining allowed region (90%CL) 
in the a — (3 plane (solid). See Fig. (7) 

Note that amongst the experimental results 
providing the SM constraints, so far the lattice 
has not entered, in any serious way, the deduc- 
tion of V ub /V cb . Of course much work is being 



done in this context too [68] but the accuracy of 
the lattice work has not yet reached the same de- 
gree of maturity as Bk or fs [69,33]. 

Figs. (7-8) show where the SM parameters 
can still lie. Hopefully improvements in the lat- 
tice calculations along with experiments will soon 
eliminate any allowed domain for the SM in these 
plots! 

Two other important quantities deducible, by 
use of the lattice results, are Vtd/Vt s and x s /xd- 

Using f B = 173 ± 40 MeV, B B = 1.0 ± .15 
along with the experimental results [40] on Bd-Bd 
mixing (xd) and the measured B-lifetime leads to: 

VtdlV ts = .22 ± .08 (43) 
Next the relation between the mixing parameters 
for B s -B s and Bd-Bd can be quantified. 



IfsJfE 



B 



B, 



x d Bb 
With the ratio f B ,/fB 
fore, R sd = B Bs /B Bd = 
eqn. (43), we get 



(44) 

= 1.16 ± .10 given be- 
1.01 ± .05, and ^ from 

1 vts 



Xs/xd ~ 18 i 14 x s /xd <C 50 (45) 
Although this ratio is far from precise the up- 
per bound appears safe and already leads to the 
important experimental implication that B s -B s 
oscillations may be observable at LEP, SLC and 
at hadron machines such as HERA-B. 

7. Summary 



1. eVe : Martinelli et al. have clarified the the- 
ory as it pertains to the lattice a good deal, 
however, relevant hadron matrix elements 
are badly needed. On some cases even 
bounds could be helpful, e.g. 06- Do these 
matrix elements suffer from large final state 
interactions? What is their a dependence? 
There is a big opportunity here for the lat- 
tice to make impact on phenomenology. 

2. B -> K*-y etc. : Data for "high m w " disfa- 
vor a constant T2 (of fo or Ai). Of course, 
T2 milder than a simple pole cannot, at 
present be ruled out. It is very impor- 
tant to quantify the q 2 dependence more 
precisely. Improved lattice calculation of 
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Rk* would yield valuable information on 
"long-distance" contributions. These would 
have important repercussions for extracting 
Vtd/Vts from measurements of B° — > p° + 7 
and B — > K*j [B^ — > + 7 is expected 
to be more "problematic"]. 

Various HQS relations work very well; not 
so V/A x . 

It may be better to use the lattice data on 
3-pt functions and 2-point functions and an- 
alyze as ratios of the two. The expected 
behavior of heavy-light form-factors (at q 2 
max) and that of the heavy-light decay con- 
stant suggests that [form factor (at q 2 max) 
/heavy-light decay constant] may have im- 
proved convergence with mht as well as 
smaller errors. 

3. Bk : JLQCD preliminary result contradicts 
Sharpe et al. result. In particular it fits very 
well to A + Ba and not at all to A + Ba 2 . 
Needs a lot of scrutiny. 

Quenching errors on Bk appear to be very 
small < 5%. All 4 studies seem to find, 

though, -Bf^ 11 to be just a bit less thanB^. 

4. Bhi Bht should be fit via Bht = const + 

Seems to give B bs (2 GeV) = .97 ± 
.06 ±.04. Again Bi, s in the full theory seems 
to be just a little less then in the quenched 
theory. Also find Bis/Ba is greater that 1 
just by a bit. 
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